Current noise and higher order fluctuations in semiconducting bilayer systems 
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We analyze the transport properties of a semiconductor based bilayer system under non- 
equilibrium conditions with special emphasis on the charge transfer statistics in the regime domi- 
nated by the exciton transport. We consider two different models. In one of them the transport 
occurs incoherently and is dominated by incoherent tunneling processes of individual excitons, while 
in the other system no disorder is present and transport processes are fully coherent. We find that 
the strength of cross correlations of currents in different layers is only insignificantly affected by 
the disorder and shows up similar behaviour in both systems. We discuss possible experimental 
realizations and make predictions for measurable quantities. 
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I. INTRODUCTION 

Electronic bilayer systems have recently drawn in- 
creasing attention due to new experimental results on 
electronic transport^— as well as due to the discovery 
of the novel macroscopic ordering driven by quantum 
effects^ - — In a generic situation such systems consist of 
two spatially separated and mutually insulated semicon- 
ducting layers which can be individually gated. On the 
other hand, the interlayer separation is engineered in 
such a way that the Coulomb interaction between ma- 
jority carriers in both subsystems is still large. This 
can be achieved by an appropriate choice for the barrier 
material. 8 ' 9 In such a situation the formation of excitons 
is possible. Even a condensation of such electron-hole 
pairs was predicted^— and subsequently observed in 
several experiments performed on GaAs quantum wells 
using AlGaAs as the insulating barrier at total filling fac- 
tor v = 1 in the quantum Hall regimeJ^— 

While the non-linear transport characteristics of such 
systems have been investigated before^, their noise prop- 
erties received only little attention. In [2l| we tried to 
close this gap and have calculated the cumulant generat- 
ing function (CGF) of a bilayer system attached to four 
different electrodes. Some very interesting effects were 
found, which can only be attributed to the non-trivial 
interactions between the layers. In order to perform the 
calculations one had to assume that the transport pro- 
cesses at different junctions occur in an incoherent fash- 
ion (sequential tunneling). Despite this severe restriction 
the currents in both layers are correlated due to the lo- 
cality of the bound electron-hole pairs. Here we would 
like to drop this assumption and consider a fully coherent 
system. 

The paper is organized as follows. In the next section 
we introduce the system and give the details of its math- 
ematical implementation. We briefly review the main re- 
sults of[2l] and adapt them to a recent experimental setup 
presented in0. In Section llHl we present a toy model con- 
sisting of two short Hubbard chains, which are coupled 
with each other, so that a formation of an interchain 



electron-hole pair or exciton is adequately described. In 
the next step we calculate the full counting statistics, 
compare the predictions of both models, analyze their 
similarities and discrepancies. We briefly summarize the 
results in Section ITVl 



II. INCOHERENT TUNNELING 
APPROXIMATION 

The basic system we study is depicted in Fig. [1] a 
semiconducting bilayer is contacted by four metallic elec- 
trodes via tunnel contacts. The Hamiltonian of the sys- 
tem reads^i 



H = H 



N ■ 



H 7 



H sc , 



(1) 



where the term describes the four metallic elec- 
trodes which are described as fermionic continua writ- 
ten in terms of electron field operators a a at chemical 
potentials [i aa . a — L,R refers to the contacts on the 
left /right side of the respective layer and a = T,B labels 
the top and bottom layer, respectively. We assume the 
leads to be in the wide flat band limit so that their den- 
sity of states (DOS) po is constant and denote the Fermi 
distribution functions in the respective electrode by n aa . 
For the phenomena we would like to consider the spin 
degree of freedom is irrelevant, therefore we work with a 
spinless system. 




FIG. 1. Sketch of the experimental setup. A semiconducting 
bilayer is contacted by four metallic leads denoted by the 
indices LT (left top), LB (left bottom), RT and RB. 
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Ht is the tunneling Hamiltonian describing hopping 
between each layer of the bilayer and the corresponding 
metallic electrodes 

H t= 7a*(a+¥„ + *+a a ), (2) 

a=T.B 

a=L,R 

using tunneling amplitudes jl.t/b, 1r,t/b- Finally, 
-ffsc describes the bilayer using a simple one-dimensional 
model as in [2l| and [22| (we shall discuss generalizations 
to higher dimensions later) 

#sc = fjx *t( x ) ^ ^ y( x ) , (3) 

where I is the longitudinal distance between the right 
and left side of the semiconductor, *S> = (^t, ^b) T is the 
two-layer spinor and Ht,Hb describe the electron/hole 
single particle Hamiltonian of the top and bottom layer. 
We describe the interlayer Coulomb interaction by an 
exciton order parameter A(sc), which in general is a 
space-dependent quantity. 8 Its absolute value Ao in 
equilibrium represents the excitonic gap and directly 
gives the exciton coherence length £sc = vf/Aq- It 
can be determined self-consistently from the interlayer 
interaction as in BCS theory^ 

In order to calculate the cumulant generating function 
(CGF) of charge transfer we use a Green's function (GF) 
method developed for quantum impurities in [24]. For ev- 
ery electrode present in the system we have to introduce 
a fictitious measuring field A = (Xlt, Alb, Xrt, Xrb) 
which counts an electron when it crosses the electrode- 
semiconductor junction in question. In this manner fluc- 
tuations, either caused by the discreteness of charge or 
thermally induced, can be taken care of. The measur- 
ing fields are both contour and time-dependent: they are 
nonzero only during the measurement time r and change 
sign on the + and — branch of the Keldysh contour. This 
procedure allows for the calculation of the CGF for ar- 
bitrary parameters in Eq. (P) as discussed in [25| and 
has already been applied to numerous quantum impurity 
problems (see e.g. 12614281) . 

The CGF can be calculated as a generalized Keldysh 
partition function 

X (X) = (T e e- i ^ TX ^ dt ), (4) 

where Tc is the Keldysh contour ordering operator and 
(t) can be derived from the tunnel Hamiltonian in Eq. 
@ using the substitution a a — > a a e lXc "*( t ^ 2 . In order to 
calculate the CGF we then rewrite the expression using 
the adiabatic potential method in the limit of long mea- 
surement times r as ln%(A) = — irU(X). The adiabatic 
potential U(X) may in turn be rewritten in terms of the 
bare and exact-in-tunneling GFs of the system described 
by Eq. ©. 



The Hamiltonian in Eq. Q is very similar to the 
BCS Hamiltonian. Indeed, T and B can be inter- 
preted as pseudospin indices to map this Hamiltonian 
with the additional exchange of B holes to electrons 
onto BCS theory^ Also one has to take into account 
that electrons and holes each have their own chem- 
ical potential use, a, & = T,B. 23 Consequently we 
can also rephrase our problem as the calculation of 
the spin-dependent CGF for a normal-superconductor- 
normal (NSN) junction with different chemical potentials 
for the spin species. We want to stress that despite the 
modelling of the condensate as a perfect homogeneous 
system in describing the tunneling we assume for the de- 
phasing time 7~0 ; £e <C VfTq <C I (incoherent tunneling 
approximation) . 

The CGFs for different normal-superconductor hybrid 
structures have been considered beforei^Sr— Here the 
task is to work out the spin-dependent case, to take care 
of the presence of electrons instead of holes in the lower 
layer and to address the case of a bilayer contacted from 
two sides. However, the A(s)-coupling in Eq. ([3]) is for- 
mally a tunneling term. Since the layers are supposed 
to be isolated from each other we have to suppress par- 
ticle current flowing between the layers. This is very 
conveniently done by adjusting the chemical potentials 
in the layers self-consistently. A similar procedure was 
implemented for a normal-superconductor-normal junc- 
tion in Ref. [35l Under such circumstances it is then 
unimportant whether we count charges on the left side 
or the right side of the respective bilayer. This means 
that we can then calculate the CGF only for the left 
side and obtain the one for the right side by exchanging 
Alt = —Xrt, All? = —Xrb- We would like to point 
out that this self-consistency condition is very different 
from the conventional one inspired by the BCS theory 
as it defines the relative position of the exciton chemical 
potential with respect to the ones in the external elec- 
trodes. 

Using the microscopic Hamiltonian approach as out- 
lined in Ref. [25|we have obtained the complete expression 
for the temperature and energy dependent CGF, valid for 
arbitrary values of the model parameters^ The expres- 
sion, which has been used for our numerical calculations, 
is quite lengthy and we do not report it here. Neverthe- 
less, all relevant ingredients of the CGF already appear 
in the limits of small bias {p-La — MHct -C A, a = T.B) 
and large bias (hlo — HRe S> A, X = T.B), where 
the expression of the CGF greatly simplifies. Indeed in 
these regimes the non-interacting self-energy due to the 
EC is either real and purely off-diagonal or imaginary 
and purely diagonal, like in a superconductor—, while 
the one due to the normal lead is always diagonal. In 
these regimes the CGF acquires the following expression 
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J2 m{l + T CT H[( ( 



a=T,B 



l)nscv(l - "ZoO] } 1 



IwJ- A 



+ ln{l + T jA («) [( 

+ (e lA ^e- lA " - l)n LB (l - n LT )] } 6 



l)n LT (l - n LB ) 

A - maxflcdrl, 



(5) 



where w t ,b = w — /isc, CT , cr 
coefBcients are given by 



T, £? and the transmission 



4r, 



(i 



r Lff ) 2 



and Ta(cj) = 



4r A 



(6) 



The energy-dependent DOS of the EC affects the hybridi- 
sations 



^La — 



TLa\u<j\ 



and = 



where T La = ir 2 Pol ct Posc7l ct A witn tne D0S of tne 
EC given by pose- n sco- refers to a Fermi distribution 
function at chemical potentials ^sc,cd <j = T,B char- 
acterising the layer of the semiconductor. We use units 
such that e = % = k B = 1 and Go = 2e 2 /h. 

The first line of the CGF describes the supra-gap con- 
tribution, which is only due to single electron transport 
in the simplified form given here and characterised by 
the normal transmission coefficient T a (oj). The second 
line describes the sub-gap contribution due to excitonic 
Andrecv reflection 22 in which an electron and a hole (in 
different layers) enter or leave the EC. psc.c: cr = T,B 
now have to be chosen such that^ 



LT 



i d 
t d\j 



XL = ((iRTh ((Ilb)) = ((IRE)). (7) 



For the sake of simplicity we consider the symmetric case 
T lt = Trx and T^b = Trb- Then the requirement 
Eq. ([7|) is always fulfilled for /zsc.T = PSC.b = if [ilt — 
-fJ*RT = -V T /2 and p L s = -(J-rb = — Vb/2. This 
condition also implies that no current is flowing between 
the upper and lower layer.— 

It is interesting to use the above formula to test fluc- 
tuation theorems. We have done so for the full result 
but for the sake of simplicity we illustrate the test for 
Eq. (J5J) - One of the most celebrated examples is the 
Cohen-Gallavotti relation, which can be directly related 
to a symmetry of the CGFj^St— For a single counting 
field A the CGF is expected to obey in the limit of long 
measurement time^ 



ln X (A)=ln X (-A-i/3y), 



(8) 



where /3 is the inverse temperature. However, in our 
case we encounter several complications compared to this 
simple case. First, the presence of an interlayer Coulomb 



I 

interaction gives rise to the subgap contribution in Eq. 
([5]) which couples the top and bottom layer and therefore 
requires to extend the simple relation in Eq. (|5]) to both 
T and B counting fields. Additionally different voltages 
can be applied on the top and bottom layer. Finally, we 
had to mind self-consistency which we satisfy by setting 
MSG.T = MSC,s = 0. Still, a relation similar to Eq. © 
holds, namely^ 2 - 



lnx(A iT , A iB ) = lnx(-A L T - i0Vr/2, -A 



LB 



iPVg/2). 



Even though this is a straightforward generalisation of 
the Cohen-Gallavotti relation in Eq. ((5|) the above re- 
lation provides an example how correlations affect the 
fluctuations in quantum systems. In this way we have 
supplied an example of a system where a generalized 
Cohen-Gallavotti fluctuation theorem holds. 

From the CGF in ((5|) we can read off several interest- 
ing transport features. At T = and \Vr\, \V B \ < 2A the 
excitonic Andreev reflections dominate the charge trans- 
port. Then we immediately obtain 



% d id 
((Ilt)} = — ttt — XL = -ttt — 
t o\ LT t a\ LB 



XL 



(I 



LB 



(9) 



Consequently, we correctly predict that the excitonic An- 
dreev reflections induce counterpropagating currents in 
the layers. 3 Depending on the bias choice we observe dif- 
ferent transport features. For equal bias in the two layers 
Vb = Vt the magnitude of the current becomes zero (ex- 
citon blockade). If only one layer is biased an opposite 
current in the other layer is induced (drag). If we further 
increase the bias difference between the layers the coun- 
terflowing currents increased This resembles the trans- 
port features observable under phase bias instead of volt- 
age bias i 22 ' 43 

So far only experiments in a typical drag configuration 
have been made. The inset of Fig. [2] (a) shows a sketch 
of the experiment described in 2\, where two Corbino 
disks are realized in a GaAs/AlGaAs double quantum 
well structure. The experiment is performed in the v = 1 
quantum Hall state of the system and we expect exciton 
formation between the two Corbino disks that serve as 
the top and bottom layer indicated in Fig. [T] In the 
top layer a current can be excited and additionally the 
induced drag current in the bottom layer is recorded. 
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We calculate the current from Eq. (|5|) and use Ft = 
0.02, T B = 0.72, A = 100 fieV. We observe good quali- 
tative agreement in Fig. [2] (a) . 

We only observe a slight deviation in the bottom current 
for voltages above the exciton gap. We ascribe this dis- 
crepancy to a small interlayer tunnel coupling that might 
still be present in the system and is also suggested from 
the analysis of the experimental data [2|. 

The next logical step is to compute the higher cumu- 
lants of the FCS. A number of details was already given in 
Ref. [U The current cross correlation of the top and bot- 
tom current for the setup shown in Fig. [5] (a) is plotted in 
Fig.[2](b) using the same parameters. As we see from Fig. 
[5] (b) one does not observe a positive cross correlation. 
This is very different from the setup in which a supercon- 
ductor is contacted to two normal drains, in which one 
observes a positive cross correlation of the two currents 
in the normal leads via crossed Andreev reflection i 33 i 44 i 45 
This discrepancy is due to the fact that we observe the 
correlation of electrons and holes and not correlations of 
electron pairs as in the case of superconductors. In this 
case the positive cross-correlation mediated by crossed 
Andreev reflection turns into an anti-correlation. We ex- 
pect that this negative correlations should be observable 
in the experiments similar to the one presented in Ref. 2 

Up to now we have considered a ID system. Due to 
pointlike contacts to bias electrodes all excitations travel 
as s-waves even in higher-dimensional systems, in partic- 
ular in genuine 2D layers. On the other hand, the ground 
state changes due to A(x) coupling between the layers are 
qualitatively the same in any dimension since the diag- 
onalization procedure is dimension-independent. In con- 
nection with the assumption of incoherent tunnelling we 
conclude that the physical picture obtained above also 
holds for 2D layered systems. We expect fundamental 
changes to occur in the case of extended contacts when 
a number of different transport channels are possible. 
However, as long as these channels are non-interacting 
the corresponding FCS would be a simple sum of expres- 
sions for individual channels. A very similar procedure 
was considered in|46| in connection with superconducting 
contacts. 



III. COHERENT TRANSPORT: A TOY MODEL 

Now we would like to abandon restrictions we made 
in the previous section. An appropriate model would 
be two non-interacting fermionic continua particle-hole 
coupled and perfectly clean so that — > oo can be safely 
assumed. Such structure is, however, quite difficult to 
handle. Therefore we model the layers by tight-binding 
chains of finite lengths. It turns out that in order to 
account for most of the transport physics we only need 
to keep two sites in each chain. The effective model can 
then schematically represented as shown in Fig. [3] 
The pairs 1,2 and 3,4 represent different layers of the 
device. The chemical potentials are now replaced by the 



bare energies of the sites ^ . Then the bilayer part of the 
Hamiltonian is given by 



H B = X] ^* fe d k d k + 7o (4 d 2 +d\d 4 ) 



k=l 



(10) 



where 70 is the tunneling amplitude between the sites 
within the respective chain and 7^ is the exciton cou- 
pling. 70 corresponds to the Fermi velocity of a con- 
tinuum realization whereas jj_ in the equivalent of the 
exciton coupling A (a;) (it is measured in different units 
though) . Coupling to the electrodes is modeled in a way 
similar to that of Eq. @ , 

H T = 7(L+c?i + L+d 3 + B^ T d 2 + R\jd 4 + H.c.) . (11) 

For simplicity we have chosen all tunneling amplitudes 
to the electrodes to be of the same magnitude and the 
DOS in each electrode to be p . In order to compute 
the FCS one performs the same transformation as shown 
in the previous sections in order to generate T x . After 
that the cumulant generating function can be evaluated 
using functional integration. We integrate the electrode 
degrees of freedom first, thereby obtaining the new ef- 
fective Green's functions for the individual sites of the 
bilayer, 



D7 



cj - fj,k - iT(2n k - 1) 



2ie lXk Tm 



~2ie- lXk T{l-n k ) -u + Mfe - *T(2n fc - 1) 



where k = 1, . . . , 4 = LT, RT, LB, RB, T = np j 2 and 
Afe are the respective counting fields. Then the CGF is 
written as 



X (A) = / V[d] e 



(12) 



where 



S = d t Kd,K 



V 



lo&z 



D 



7J-CT 2 




\ 

7-L cr z 

10° z 

D-, 1 j 



(13) 



d are superfields defined by 

d = (d 1 - 1 di +1 d 2 -,d2+,d 3 - 



,d 4+ ). (14) 



Here the index ± denotes the Keldysh component of the 
respective field. After the calculation of the functional 
integral we finally obtain 

lnx(A) = TG J rfwln[detK(A)/detK(A = 0)] (15) 

lik are yet to be appropriately chosen. To fix them we 
perform the same procedure as in Section ITT] We com- 
pute the currents on both junctions of the upper layer 
and equalize them. On the other hand, from symmetry 
arguments it turns out that the choice /it = Ms = is 
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100 200 300 
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FIG. 2. (a): comparison of the prediction for the current using Eq. ([5]) for a typical drag experiment in ECs. We use fitting 
parameters Tt = 0.02, Tb = 0.72, A = 100 fieV and take the temperature T = 17mK from experiment. The dashed blue and 
red curve correspond to the top and bottom current data from the experiment, respectively and the dashed black curves are 
the results from our model. 

(b): noise for the top current (upper blue curve) and the cross correlation of the top and bottom current (lower red curve) 
using Eq. ((5]) and the same fitting parameters as in (a) for the experiment 0. 



Yo 



Yo 



FIG. 3. Sketch of the toy model for coherent transport. The 
semiconducting bilayer in Fig. [1] has been replaced by two 
pairs of tight binding chains which have an excitonic coupling 
described by -y± . 



particularly convenient. One interesting parameter con- 
stellation is the situation of no bias in one of the layers 
(the upper one for definiteness) and finite voltage applied 
symmetrically around fis in the bottom subsystem. One 
can produce analytical formulas for the direct current in 
the bottom layer (v = d) as well as for the induced cur- 
rent in the top layer (y = ind), valid at all temperatures, 



G f 

— / du)D v (w)[n LB (u) - n RB {io)] , 



(16) 



where n a (uj) are the Fermi distribution functions in the 
respective electrodes. The effective transmission coeffi- 
cients Dv(uj) have an interesting structure. They are 
particularly simple under the assumption 70 = 71 =7/2, 



D d (u) 



And(w) 



7 2 r 2 [ 7 2 + 3(r 2 



■a, 2 )] 



(r 2 + w 2 )[ 7 4 + 2 7 2 (r 2 - oj 2 ) + (r 2 + u 2 ) 2 } 

7 2 r 2 [ 7 2 -(r 2 + ^ 2 )] 

(r 2 + w 2 )[ 7 4 + 2 7 2 (r 2 - w 2 ) + (r 2 + w 2 ) 2 ] 



First of all we observe that Dd(uj) is a positively defi- 
nite quantity and has no zeros, as one would expect from 
a conventional transmission coefficient. Being a 'good' 
transmission coefficient, its denominator does not have 
any zeros on the real axis either. The induced current 
is finite in almost all parameter regimes. However, for 
7 2 > T 2 And(^) changes sign at w = y/j 2 — T 2 . As a 
result two different generic behaviour types can be ob- 
served. For weak I7I <C T imd(V) flows in the direction 
opposite to Id(V). This regime is largely compatible with 
the exciton blockade case and we shall call it 'exciton' 
regime. Here the hybridization of the layer states with 
their 'own' electrodes is dominant. In the opposite case of 
large 7 2 > T 2 + V 2 the currents flow in the same direction 
and one is confronted with the standard drag situation. 
These two regimes can also be observed also for 70 =/= 7 j_ , 
see Fig. HJa). Strong 7 couplings tend to destroy the ex- 
citonic state favoring a rearrangement of individual en- 
ergy levels in such a way that from the point of view of 
incoming electrons the constriction is a larger quantum 
dot rather than two individual layers. Nonetheless, the 
particle countcrflow might be significant. This is the fea- 
ture our toy model shares with the system considered in 
the previous Section [TT1 

In the next step we compute the cross correlation of 
currents in the layers. The calculations are as straight- 
forward as in the case of the current. The results for 
the cross cumulant of the induced and direct currents for 
representative parameter sets are shown in Fig. H](b). It 
turns out that the correlations are remarkably strong and 
as far as their amplitude is concerned their voltage be- 
havior quantitatively follows the results of the incoherent 
approximation. A fundamental difference is however, the 
sign of the correlations which never changes sign and is 
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FIG. 4. Plots of the current and the cross correlation of currents, (a) shows the current 7 as a function of voltage V for 70 = F 
and T = 0.05r. The solid black curve corresponds to the induced current for an excitonic coupling y± = 0.85F while the direct 
current is shown as the dotted dashed curve. The dashed black curve refers to the induced current for 7x = 4r and the direct 
curve is shown as the dotted curve. 

In (b) we show the cross correlation of the direct and induced current for the same temperature and 70. The solid black curve 
is for 7x = 0.85r, the dotted dashed black curve shows the result for y± — 2T, the dotted curve refers to 7^ = 3r and the 
dashed black curve represents y± = 4T. 



always positive, unlike the situation discussed in Section 

nn 



IV. SUMMARY 

As we have seen above, a drag counterflow due to the 
exciton coupling is present in both models. This is a very 
remarkable result in view of the fact, that in the first 
model transport is assumed to be completely incoherent 
(t^ is finite and small). On the other hand, while the 
cross cumulant of the induced and direct currents turns 
out to be positive in the second model, it is negative in 
the incoherent approximation. From the physical point 



of view the fundamental difference between the models 
is the fully developed gap in the excitation spectrum of 
the incoherent model. It cannot be generated in the toy 
model though. Therefore we conclude that for our setups: 
(i) for the observation of the drag counterflow the gap is 
not necessary while it is necessary for perfect Coulomb 
drag; (ii) in the presence of the exciton coupling the coun- 
terflow can be achieved even in disordered systems; (iii) 
current anticorrelations exist only in the gapped systems. 
We expect these predictions to be observable in the up- 
coming experiments. 
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F. Dolcini for many interesting discussions. The finan- 
cial support was provided by CQD of the University of 
Heidelberg. 



R. V. Gorbache v et al., Nat. Phys. 8, 896 (2012). 

D. Nandi et al., |Nature 488, 481 (2012)| 

J. Su and A. H. MacDonald, Nat. Phys. 4, 799 (2008). 

J. P. Eisenstein and A. H. MacDonald, Nature 432, 691 

(2004). 

A. A. High e t al., N ature 483, 58 4 (2012). 
Z. W. et al., |New J. Phys. 14, 063010 (2012)) 
A. Perali et al., ArXiv e- prints 



arXiv:1211.6562 [cond-mat.supr-con] 



(2012), 



Phys. Rev. B 78, 121401 (2008) 



I. 



Shevchenko, 



H. Min et al.. 

A. I. Bezuglyj and 
|J. Low Temp. Phys. 37, 583| 
1 D. Snoke, Science 298, 1368 (2002). 
Y. E. Lozovik and V. Yudson, Zh. ksp. Teor. Fiz. 71, 738 
(1976). 

J. M. Blatt et al.. 



A. 



A. 



Phys. Rev. 126, 1691 (1 962) 



Pikalov 



and 



D. 



V. 



Fil, 



|Nanoscale Res. Lett. 7, 145 (2012)] 

J. P. Eisenstein e t al., |Phys. Rev. Lett. 68, 1383 (1992) 
M. P. Lilly et al., |Phys. Rev. Lett. 80, 1714 (1998jT~ 
I. B. Spielman et al., |Phys. Rev. Lett. 84, 5808 (200"0)| 



M. Kellogg et al., Phys. Rev. Lett. 93, 036801 (2004) 

E. Tutuc et al., |Phys. Rev. Lett. 9 3, 036802 (2004JP 
L. Tiemann et al., New J. Phys. 10, 045018 (2008) 

A. D. K. Fink e t al.,|Phys. Rev. Lett. 106, 236807 (2011) | 
H. Soller et al., |Phys. Rev. Lett. 10 8, 156401 (2012) 

F. Dolcini et al., |Phys. Rev. Lett. 104, 027004 (2010)| 
C. 



Comte 



and 



|J. Phys. France 43, 1069 (1982) 



O. 



Gogolin 



and 



A. 



Nozieres, 



Komnik, 



Phys. Rev. B 73, 195301 (2006) 



('. Cu.-vhs <-r al.. Phys. Rev. L! 54. 7:5(>(> i 1996) 



J 

B. Dong et al., |Phys. Rev. B 84 
R. Avriller and AT 
|Phys. Rev. B 80, 041309 (2009) | 
D. Breyel and 

|Phys. Rev. B 84, 155305 (2011)| 
Y. E. Lozovik and O. L. Berman, Zh 
1879 (1997). 
' W. Belzig and Y. 
|Phys. Rev. Lett. 87, 067006 (2001)| 



235408 (2011) 
Levy 

A. 



Yeyati, 



Komnik, 



ksp. Teor. Fiz. Ill, 



V. 



Nazarov, 



7 



B. A. Muzykantskii and D. E. Khmelnitskii, 
|Phys. Rev. B 50, 3982 (1994) | 
D. Futterer et al., |Phys. Rev. B 79, 
H. Soller and A. Komnik, Eur. Phys. J. 
H. Soller and A. Komnik, |Physica E 44, 425 (2011)] 

C. J. Lambert, J. Phys.: Condens. Matter 3, 6579 (1991). 
R. Fazio and R. Raimondi, 



05450 5 (2009)] 
D 63, 3 (2011). 



Phys. Rev. Lett. 8 0, 2913 (1 998) 



S. 



Hershfield 



and 



H. 



Zhao, 



Phys. Rev. B 56, 3296 (1997) 



D. J. Evans et al., |Phys. Rev. Lett. 71, 2401 (1 993) 



G. Gallavotti and E. G. 

|Phys. Rev. Lett. 74, 2694 (1995) | 

A^ Komnik and H. 



41 

42 



D. 



Cohen, 



Saleur, 



Phys. Rev. Lett. 10 7, 100601 (20 11) 
M. Esposito et al., |Rev. Mod. Phys. 81, 1665 (2009) | 
K. Saito and Y. Utsumi, |Phys. Rev. B787 115429 (2008) 
S. Peotta et al., |Phys. Rev. B 84, 184528 (2011)1~ 
P. Recher et al., |Phys. Rev. B 63, 165314 (2001)] 
J. Wei and V. Chandra sekhar, Nat. Phys. 6, 494 (2010). 



